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Abstract 

We examine discrete gauge symmetries in axionic extensions of the SSM which 
provide a solution of the //-problem. Automatic-PQ symmetry and proton sta- 
bility are shown to be guaranteed by certain discrete symmetries. Focusing on 
the L-violating discrete symmetries we discuss two sources of neutrino masses 
and their relevance for the solar neutrino problem. 
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In the conventional standard model baryon (B) and lepton number (L) invariance 
are automatic consequences of the gauge symmetry and the particle contents. Contrary 
to this the supersymmetric standard model (SSM) has to be endowed with certain 
discrete symmetries (typically R parity) to enforce the B- and/or L-invariance. An 
interesting possibility is to interpret these discrete symmetries as gauge symmetries 
which can be remnants of a broken unifying gauge group. Moreover, discrete gauge 
symmetries can stabilize B- and L-invariance even if quantum gravity effects violate 
global symmetries. The discrete symmetries should then fulfill certain anomaly-free 
conditions j^, ^. Recently, it has been found that proton stability or solar neutrino 
mixing can be consistent with certain classes of anomaly-free discrete symmetries 0, ^ . 

In this paper, we address the question of an "automatic" PQ-symmetry in the 
SSM guaranteed by a certain discrete gauge symmetry which is assigned to prevent 
serious B- or L-violation. As a consequence all desirable global symmetries of the 
model can be made stable against possible quantum gravity effects. 

As we will see it appears to be very difficult to implement an automatic PQ- 
symmetry of the Kim- type 0. Therefore we concentrate on the Dine-Fishler-Srednicki- 
axion ||^ and we extend the Higgs sector of the SSM by adding some standard model 
singlets Ai (PQ-fields). The requirements of automatic PQ lead to a solution of the 
/i-problem Q by introducing the d = 5 operator HHA^. When the field A gets a 
vacuum expectation value (VEV) of the order of {fa is the PQ-symmetry breaking 
scale with 10^'' GeV < fa < 10^^ GeV) the generated /i-term receives the correct order 
of magnitude. 

Depending on the PQ-field contents and discrete symmetries, various kinds of B-, 
L- or PQ-violating operators appear. One should check whether some combinations 
of these operators can be eliminated to ensure proton stability and the axion solution 
of the strong CP problem. In particular, the axion solution is very sensitive to the 



presence of higher dimensional PQ-nonsinglets because of the small axion mass |TD 



The fermionic part of the chiral superfields Ai can serve as a right handed neutrino 
(RHN) due to the possible couplings LHA and LHA^. Clearly neutrino masses can be 



generated by nontrivial VEV's < A>. In addition, the see-saw mechanism [|TT| can be 
operative depending on the masses of the fields A. Indeed this amounts to the super- 
symmetric version of the neutrino-axion connection investigated in \V2\. Another pur- 
pose of this paper is to discuss the phenomenological implications for neutrino physics in 
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this context. Especially we address to possible solutions of the solar- neutrino-problem 
recently confirmed as a problem of elementary particle physics by the measurement of 
the GALLEX-Collaboration [0, g. 



We begin our discussion with the DFS-type axionic extensions of the SSM. 
The classification of discrete gauge symmetries for the SSM was given by Ibanez and 
Ross [0]. They wrote all discrete symmetries qn compatible with "standard" operators 
in terms of three generators : 

= R^A^L' (1) 
which leads to charges of the standard particles as given in table |1]. Additionally this 
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Table 1: Discrete charges of standard particles. 

symmetries have to fulfill anomaly cancelation conditions due to the gravity-, SU (3)- 
and S'f/(2)-interaction. In our extension the generators of eq. (|ip have to be continued 
to the new fields Ai. Moreover new generators can appear which only act non trivially 
on the Ai. What happens in detail of course depends on the explicit form of the 
superpotential. Some examples will be given later. The modified form of the anomaly 
cancelation conditions reads 

3m + n + k = -rjsN for gravity 

3n = for5f/(3) (2) 

31 + n = for5?7(2), 

where 77 = 1 for even and r] = for A^ odd. Since we added SM-singlets only the 
gravity anomaly is modified hj k = J2(^{Ai)- Parts of our discussion can be done 
referring to the value of k instead of using the explicit form of W{Ai). 

We start with the implications of the automatic-PQ idea. An explicit HH-term 
violates the PQ-symmetry and so we should demand n 7^ 0. The last two anomaly free 
conditions show that only a restricted class of discrete symmetries is compatible with 
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this requirement 



N = 9r, n = 3r, I = 2r, r = 1,2,3,... (3) 

Having forbidden HH we are forced to generate this term dynamically. This amounts 
to solving the /i-problem. Introducing the conventional DFS-coupling i^^fyli, however, 
needs a very small coupling constant. Therefore we use one of the terms HHA1A2 or 
HHA\ which naturally lead to the desired order of magnitude. For discrete symmetries 
with = 3* the second term cannot be used because the coupling A\ then violates 
PQ-symmetry. Moreover PQ-symmetry has to be respected by higher dimensional 
operators. As shown recently [l^, the lowest PQ-nonsinglet operators constructed out 



of only PQ-fields should be of (component field) dimension d = 10. But the restriction 
is less stringent for the standard fields since they have much smaller VEV's than PQ- 
fields. For standard operators, it is enough to forbid LH [], i. e. m 7^ r, to protect the 
axion solution. For mixed operators (operators containing standard fields and fields 
Ai) only the combination of LHAi and a PQ-nonsinglet LHAjAk can cause a problem. 
Of course this depends on the VEV's of the A's. As we will see from the examples, it 
is possible to forbid either of these two if we do not want to introduce too many fields. 
Dangerous d = 7 operators generally appear in the pure Aj-sector. To make a concrete 
discussion, let us take a /i-term HHA\ and a superpotential W = {A1A2 — fa)Y. As 
a consequence of the second anomaly-free condition in eq. (0), one finds the allowed 
operators : A^ or etc. Now we have to worry about the PQ-nonsinglet operators 

in the scalar potential like afog?/*. It can be obtained from the terms A1A2Y and Af 
in the superpotential. At first glance, it does not make a problem since < Y >= in 
the global SUSY case. But if we take account of the full supergravity Lagrangian, the 
field Y can develop a VEV of the order of f^^^/M^i. Here p can depend on the specific 
form of the superpotential |1^. For instance, the above superpotential gives p = 1. 



So, one has to take a superpotential which gives p > 3 to protect the axion solution. 
From this, we favor a very light axino (mass < 10~^eV) since its mass is proportional 
to < F >. 

Now we turn to B-, L-number. The discrete gauge symmetries can be classified 
according to the allowed pattern of the B-, L- violating standard terms LQd, LLe and 



^This is also necessary in the spirit of the solution to the //-problem. 
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forbidden 


allowed for some m 


fc- values for = 18 


c = mod 3 


L- violating 


B- violating 


3,-6 


c = r mod 3 


B- violating 


L- violating 


6,-3 


c = 2r mod 3 


all 


none 


0,9 



Table 2: Patterns of allowed and forbidden operators for different classes of potentials 
W{A,). 



udd 



GMP : m ^ 4r, 6r 
GLP : m = 6r 
GBP : m = 4r . 



(4) 



Proton decay safety on higher dimensional standard operators is ensured if QQQL, uude 
and the combinations {udd, QueH/QuL* /ud*e) and {LQd, QQQH) are forbidden. This 
is fulfilled for all our symmetries in eq. ^ ! Additionally new B-, L-violating opera- 
tors can appear by adding singlet fields Ai. The combinations {LQdA^, uddA'j) with 
p + q < 3 should be nonsinglets under discrete symmetry to avoid fast proton decay. 
Interesting patterns of allowed and forbidden B- and L-violating operators with dimen- 
sion < 5 appear if all charges a{Ai) are multiples of 3. This is always the case if the 
coupling HHA\ is demanded to be present. But also for HHA1A2 there exists a class 
of superpotentials with this property. We define c = m mod 3. Then for r mod 3 7^ 
three patterns of operators occur two of them automatically lead to a proton decay 
safe theory (see table 0). For r mod 3 = all operators can be allowed in the case 
c = whereas in the opposite case all operators are forbidden. The /c-values in terms 
of the class c are given by 

k 



mod 3 



r even 



(5) 



3 I — r — c r odd 

We saw that certain discrete symmetries can allow L-violation. In this class of 
symmetries neutrino masses can be generated with the fields Ai serving as RHN's. 
The relevant couplings are the following : 



anLiHA. + ^LiHA.Ak + ^HHAiA^+ 



4 



LjHLiH H LjH H H . 

2M ^ M 



The last two terms - if present - can produce neutrino masses O(10~^) eV (the last 
term only one mass of this order). With this in mind we neglect them in the following 
diagonalization procedure. The mass matrix of / = {ui, h,h,dk) can be read off from 
the above couplings : 

/ c Md \ 
fi 
^ /i 
\Ml h a Mm ) 

Neglecting again contributions O(10~^) eV one can diagonalize this matrix in three 
steps : rotating away the v — /i-mixing, diagonalizing Mm and applying see-saw mech- 
anism for all "large" eigenvalues of Mm- Then we arrive at : 
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(6) 
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The diagonalized matrix Mm has been divided into three parts : the axion and very 
small eigenvalue (<C ^3/2) psnt A, the small eigenvalue part 5* with eigenvalues ~ ^3/2 
and the large eigenvalue part L with eigenvalues 0{fa)- The other matrices have been 
divided into the corresponding parts. The equations show that the couplings LiHAj 
may only be present if the VEV of Aj is small (~ 0(7713/2)) and Aj does not mix with 
the light eigenmodes in the S- and A-sector. Under these assumptions the matrices 
Da and Ds are at most of the order /am3/2/M ~ keV. The application of the see-saw 
mechanism for Ds is then justified. From the matrix (0) the following conclusions can 
be drawn : 



1. In general Da will be nonzero if couplings LHAA are present. Then the axino 
mass determines the mass of one neutrino mode. For a light axino mass (<C 
keV) a keV-Dirac-neutrino (z/, a) appears. Though there are some hints for a 17- 
keV-neutrino experimental results as well as cosmological implications 
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disfavor this possibility. For an axino mass in the keV-range a keV-Majorana- 
neutrino is present which is excluded by neutrino-less double beta decay. It 
remains the possibility of a heavy axino (0(10^ GeV)) generating a Majorana 
neutrino mass of O(10~^) eV. This is the correct order for the MSW-solution 
of the solar neutrino problem. However, the last possibility which is favorable 
from the viewpoint of neutrino physics is in danger to spoil the axion solution as 
discussed. 

2. Couplings LHA can be present under the above requirements. Together with the 
see-saw mechanism they lead to 0(eV) Majorana neutrinos (for couplings 0(1)). 

3. Apart from these mechanisms neutrino masses relevant for the solar neutrino 
problem can be generated only by see-saw suppression through light fields Ai : 
If Ds is of 0{fam3/2/M) ~ keV then Ai masses of 0(7713/2) lead to masses 
0{flm^/2/M^) ~ 10-3 eV (for large fa ~ 10^^ geV). A possible source 
of a mass of 0(7/13/2) can be radiative corrections |TP| . 

4. The couplings LHLH cause Majorana neutrino masses of 0(10^^ eV) in general 
for all flavors. The coupling LHHH can generate only one mass of this order. 

5. Neutrinos with mass 0{eV) which can be candidates for Dark Matter can be 
achieved in two ways. According to the first remark the scale fa ~ 10^^ GeV 
can combine with an axino mass (~ 1 GeV). Alternatively the superpotential can 
allow for a see-saw suppression by a scale of 0{fa) as in the second remark. 

Remarkably solving the solar-neutrino-problem via masses as discussed or via LQd j^/ 
is only possible for the second pattern in table Proton decay safety follows then 
automatically. 

Before going to the examples let us discuss the relation of the automatic PQ- 
idea and the Kim-type realization of axion. The essential coupling for the Kim-axion 
K1A1K2 (where Ki and K2 are heavy quarks) leads to the following parameterization 
of charges : Ci{Ki) = —p -\- q, a{K2) = —p — q, a{Ai) = 2p. The anomaly cancelation 
conditions read 



3m + n + 2p + k 



3n -|- 2p 
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A difficulty in implementing the Kim-axion is that the model is in danger to possess 
an additional DFS-PQ-symmetry. Then the introduction of additional quarks and 
moreover the whole Kim-PQ-mechanism is superfluous. One possibility to avoid this 
is to forget about solving /i-problem and demand HH to be present. But in this case 
n = and the anomaly free conditions imply a{Ai) = 0. So there is no chance for 
automatic PQ. On the other hand, if 7^ at least one coupling T which is invariant 
under discrete n-symmetry but noninvariant under the corresponding global U{1) is 
required : a{T) = sn = s'N (s 7^ 0) and a[/(i)(T) 7^ 0. Now for all small values of 
s {s = 3/2,2,5/2..) the anomaly cancelation conditions show that some small power 
of Ai is allowed. But the generation of large s-values leads to complicated potentials 
and the introduction of many fields Ai. We conclude that an automatic pure Kim-PQ 
extension of the SSM cannot be implemented in an easy way. 

Now we discuss all acceptable potentials W for N = 9, 18 up to four fields A^. 
Later we will give one example for = 36 which naturally can allow three couplings 
LHNi 0. The superpotential W should have the following general properties : The 
superpotential is the most general one compatible with its discrete symmetry. Of 
course it should admit a PQ-symmetry. The discrete charges a{Ai) should satisfy the 
constraint = k. To force the PQ-fields responsible for the //-term to get a 

nonvanishing VEV the existence of a field Y with a{Y) = is necessary. Finally, 
supersymmetry is supposed to be unbroken in the sector of singlet fields A^. One can 
assume the Polony form for the supersymmetry breaking. 

Let us ffist consider the case with = 9. Here we have r = 1. We allow the coupling 
HHA1A2 assuming the discrete charges of the A's to be a{Ai, A2) = (—1, —2) 0. The 
anomaly free condition in eq. (|^) allows the solutions k = 3x for an integer x given by 
m = 2 — X (mod 3). We find that our minimal choice for the superpotential requires 
two more fields with 0(^3, A4) = (1, 2). Dropping the cubic and quadratic term for Y, 
the corresponding superpotential reads : 

W = y.AlA, + y^A^Al + (y^A^A, + y^A^A^ - f)Y (9) 

where f ^ fa- This corresponds to A; = 0. 

^AU PQ-fields with this kind of coupling will be denoted by N in the following. 

■^Two other discrete charge assignments are possible : a{Ai, A2) = (2, 4) or (—4, 1). The following 

arguments are also applicable to these cases. 
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From the previous discussion we know that certain combinations of B- and L- 
violating operators should be avoided {{LQdA^ ,uddA'^) with p + g < 3, LH , LHA). 
For any m at least one of these operators is allowed. Since we can hardly expect 
that introducing more fields can avoid this difficulty, we exclude the possibility for Zg 
symmetry. 

The next examples are the Zig symmetries with r = 2. The values of k,m follow 
the pattern k = 3x, m = 1 — x (mod 3) for an integer x. We now introduce a PQ-field 
Ai with discrete charge —3 to allow the coupling HHAf. Automatic PQ-symmetry 
requires a{Ai) ^ ±6 for all PQ-fields. In agreement with table 0, we find three classes 
of symmetries which share the similar phenomenological property : 

(I) m = 1 mod 3; /c = 0, 9 : The minimal choice of the PQ-fields can be achieved by 
taking two fields Ai,A2 with a{Ai,A2) = (—3,3). Then we take the superpotential 
W = {\A1A2 — P)Y + itinN'^, where f ^ rriN ^ fa- The presence of the field N is 
necessary to achieve /c = 9. One finds that this class of discrete symmetries allows no 
d = 5 operators. Therefore, B- and L-symmetry are insured quite strongly. 

(II) m = mod 3; A; = 3, —6 : For this case we need three fields with a{Ai, A2, A3) = 
(—3,3,3) for k = —3 and one more field for k = —6. The corresponding superpo- 
tential can be read off from eq. ( pUj ) by exchanging Ai and A3. As expected, this class 
of discrete symmetries respects L-invariance. But B-violation can appear in renormal- 
izable terms or in the higher dimensional operators with PQ-fields. For instance, the 
discrete symmetry allows udd, uddN, Q^H, uddA2^3 and uddAi for m = 8, 3, 6, 9 and 
—3 respectively. 

(III) m = 2 mod 3; k = —3, 6 : In this case, we have three fields with a{Ai, A2, A^) = 
(-3,-3,3) for k = —3 and one more field A^ for k = 6. One can build two superpo- 
tentials satisfying our general requirement, 

\ (A1A1A3 - /2)Fi + (A2A2A3 - /|)r2 + m^vA^' 

The first superpotential has two massless multiplets : \2Ai — A1A2 and the axion 
multiplet. In the global SUSY case, the VEV of the first multiplet is not determined. 
But it is set to be if a supergravity Lagrangian with minimal kinetic terms and a 
vanishing cosmological constant is considered. The second one gives nontrivial VEV's 
to all the PQ fields and only one massless (axion) multiplet is present. These features 
become important for neutrino physics. 
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Contrary to the above case, B-invariance is now insured. We discard the cases with 
m = 2 and 5,-1 since they allow either LH or LHAi. The L- violating operators are 
LQd, LLe and LHAl ^ for m = 8; LHAl and LHHH for m = -4; LHN, LQdAi^2, 
LLeAi 2 and LHLH for m = —7. According to the general discussion the neutrino 
masses in this three cases behave as follows : 

(i) m = 8, —4 : As discussed before, the field Y can get a VEV of order of 10^ GeV to 
produce a heavy axino. Then, there may appear one neutrino with mass O(10~^) eV for 
the both cases in eq. (|10|)) ^"^^ one keV-Dirac neutrino for the first case only. However, 
this large < Y > spoils the axion solution. To avoid this, we need a superpotential 
with sufficiently small gravity-corrected VEV < Y >. This may result in two (one) 
keV Dirac neutrino for the first (second) superpotential. In this minimal example, we 
cannot, therefore, achieve the desirable features for neutrino masses. 

(ii) m = —7 : If the field is present {k = 6) it is the only one which can couple to 
LH. It leads to one Majorana neutrino with mass 0(eV). Here we can take a small 
VEV < y > to protect the axion solution without arriving at a keV-Dirac neutrino. 
Due to the allowed coupling LHLH the other neutrino modes will have in general mass 
0(10-5) eV. 

Up to now, we discussed the simple case with one PQ-field coupled to HH. If the 
term HHA1A2 is forced to be allowed we find that - as in the Zg-case - at least four 
PQ-fields are necessary to build the potential. For a{Ai) = 3, 9, 15 all charges a{Ai) 
are multiples of three and the classification of operator patterns from the general part 
is applicable. We find four types of minimal potentials which cover the three possible 
patterns of allowed and forbidden operators. For a{Ai) ^ 3, 9, 15 the only minimal 
potential is given by eq. (P) with Zig-symmetry. Because of the extended symmetry it 
cannot be ruled out as in the Zg-case. 

We saw that the most promising way of generating neutrino masses without spoil- 
ing the axion solution is to introduce couplings LHN . In eq. (|TD|) the mass of the 
right-handed neutrino was put in by hand. Now we construct an example which 
shows a relation between the axion scale and the right-handed neutrino scale. For 
this, let us go to the Z36 symmetry. We take PQ-fields whose discrete charges are 
given by a(y4i^2, ^3, A'1^2,3) = (—6,6,3). This allows us to build a superpotential 
W = (XiAiAs + X2A2A3 - p)Y + h.A^Nl with k = Now the discrete symme- 
try with m = 1 permits the desired coupling LiHNj and forbids any d = 5 operators 
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as well as B- and L-violating d = 4 operators. One finds that the masses of right- 
handed neutrinos are given hj hi < >■ Irrespectively of < F > the axino receives 
a radiative-corrected mass of O(IO^) GeV by the coupling A^N^ . 



We finish our discussion with some remarks about generation dependent discrete 
symmetries. The leptonic charge corresponding to / is split into three generation parts 
li. Only the SU{2) anomaly cancelation condition changes : L + n = with L = Y^li. 
The "automatic" requirement n leads to = 3r, n = r and L = —r. So gener- 
ation dependency lowers the smallest N which admits automatic PQ to = 6. The 
potentials for = 6 are our = 18 potentials from the generation independent case. 
As in the generation independent case HHA\ is forbidden for A^ = 3*. Since n is not 
necessarily a multiple of 3 any more the operator classification for the term HHA\ is 
not valid. Only if A^ = 9r this is still satisfied and general conclusions can be made : 
There exist two types of (mod 3)-combinations of the U : 
generation independent : (/i, I2, h) mod 3 = (0, 0, 0), (1, 1, 1), (2, 2, 2) 
generation dependent :(Zi, I2, 13) mod 3 = (1, 2, 0)+perm. 

Apart from Q^L and u'^de which were forbidden in the generation independent case 
table 1^ remains valid for the B- violating operators. The L-violating operators with only 
one leptonic family index are allowed for m+li mod 3 = 0. Eq. gives the correspond- 
ing values of k. Consequently for the second type only one neutrino flavor can receive 
"large" masses by one of the mechanisms discussed. Interestingly this generation de- 
pendent symmetries differ from the so called type II symmetries (i. e. h = h h) 0] 
necessary to solve the solar neutrino problem by allowing Li ^^Qd. 

In conclusion we discussed axionic extensions of the SSM with a discrete gauge 
symmetry to stabilize PQ-, B- and L-invariance. Our requirement of an automatic PQ- 
symmetry excluded the Kim-type axion. This led us to implementing the DFS-type 
axion and to solving the /i-problem. We showed that a certain class of Zj^ symmetries 
with A^ = 9r can ensure proton stability and simultaneously protect the axion solution. 
To suppress the effect from higher dimensional operators a small tree-level axino mass 
(< O(10~^) eV) is favored. For a yU-term HHA\ we found three classes of discrete 
symmetries which show characteristic patterns of B- and L-violating operators. One 
pattern of them allows only L-violating operators 
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For this class two sources of neutrino masses, namely LHAA and LHN, were 
discussed. They represent different types of neutrino-axion connections with RHN 
scales m3/2 and fa, respectively. If LHAA is present in general one neutrino mass is 
fixed by the mass of the axino. This excludes an axino mass O(keV) resulting in a 
keV-Majorana neutrino. An axino mass of 0(m3/2) can generate a favorable neutrino 
mass of O(10~^) eV though the axino solution is in danger. On the other hand a low 
axino mass leads to a keV-Dirac neutrino. Therefore we preferred the coupling LHN. 
A term AN"^ in the superpotential then can produce a RHN mass rriN fa leading to 
Majorana neutrinos with mass of 0(eV). The axino receives a one loop radiative mass 
~ 777,3/2 which contrary to the tree level mass does not effect the axion solution. One 
model with this aspect was constructed for Z^q- 

For the two simplest symmetries and Zi^ we discussed all acceptable superpoten- 
tials with less than five PQ- fields A^. The symmetry turned out to be incompatible 
with our requirements. In the Z\^ case we gave examples for each of the three operator 
patterns and the above aspects of neutrino-axion connections. 
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